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1 Introduction and preliminaries 



1.1 Throughout this paper, we shall assume that X is a linear regular re- 
current diffusion taking values in R + with an instantaneously reflecting 
boundary. Let and E x denote, respectively, the probability measure and 
the expectation associated with X when started from x > 0. We assume that 
X is defined in the canonical space C of continuous functions u> : R+ i— ► R + . 
Let 

Ct := cr{u)(s) : s < t} 

denote the smallest cr-algebra making the co-ordinate mappings up to time t 
measurable and take C to be the smallest cr-algebra including all cr-algebras 
C t , t > 0. 

The excursion space for excursions from to associated with X is a 
subset of C, denoted by E, and given by 

E := {e £ C : e(0) = 0, 3 £(e) > such that e(t) > V t G (0, ((e)) 

and e(t) =0Vt> £(e)}. 

The notation Et is used for the trace of Ct on E. 

As indicated in the title of the paper our aim is to gather a number of 
fundamental results concerning the excursion theory for the diffusion X. In 
Section 2 the classical descriptions, the first one due to ltd and McKean and 
the second one due to Williams, are presented. In Section 3 the stationary 
excursions are discussed and, in particular, the description due to Bismut 
is reviewed. After this, in Section 4, we proceed by analyzing excursions 
straddling an exponential time. The paper is concluded with an example on 
Ornstein-Uhlenbeck processes. 

Our motivation for this work arose from different origins: 

• First, we would like to contribute to Professor Ito's being awarded the 
1st Gauss prize, by offering some discussion and illustration of K Ito's 
excursion theory, see [21j . when specialized to linear diffusions. The 
present paper also illustrates Pitman and Yor's discussion (see [42] 
in this volume) of K. Ito's general theory of excursions for a Markov 
process. 

• In the literature there seems to be lacking a detailed discussion on 
the excursion theory of linear diffusions. Information available has a 
very scattered character, see, e.g., Williams [52], Walsh [51], Pitman 
and Yor [38], [39], 00], 01], Rogers 06], Salminen 09]. The general 
theory of excursions has been developed in Ito [21], Meyer [35], Getoor 
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03], Getoor and Sharpe [E], [H], [E], [T7J, Blumenthal Although 
the case with Brownian motion is well studied and understood, for 
textbook treatments see, e.g., Revuz and Yor [43J and Rogers and 
Williams [17], we find it important to highlight the main formulas for 
more general diffusions using the traditional Fellerian terminology and 
language. 

• To generalize some recent results (see Winkel |53] and Bertoin, Fujita, 
Roynette and Yor [2]) on infinite divisibility of the distribution of 
the length of the excursion of a diffusion straddling an independent 
exponential time. 

• The Ornstein-Uhlenbeck process is one of the most essential diffusions. 
To present in detail formulae for its excursions is important per se. 
One of the key tools hereby is the distribution of the first hitting time 
H y of the point y from which the excursions are observed. For y = 
this distribution can be derived via Doob's transform (see Doob [8j) 
which connects the Ornstein-Uhlenbeck process with standard Brow- 
nian motion (see Sato [50], and Going- Jaeschke and Yor pj5]). For 
arbitary y the distribution is very complicated; for explicit expres- 
sions via series expansions, see Ricciardi and Sato [H], Linetsky [32] 
and Alili, Patie and Pedersen [T]. We will focus on excursions from 
to and relate our work to earlier papers by Hawkes and Truman [20J , 
Pitman and Yor [40] . and Salminen [IB]. Due to the symmetry of the 
Ornstein-Uhlenbeck process around 0, it is sufficient for our purposes 
to consider only positive excursions - the treatment of negative ones 
is similar - and view the process with values in R + and being a 
reflecting boundary. 

1.2 In this subsection we introduce the basic notation and facts concerning 
linear diffusions needed in the sequel. A main source of information remains 
ltd and McKean [22J, see also Rogers and Williams [17], and Borodin and 
Salminen [6]. 

(i) Speed measure m associated with X is a measure on R + which satisfies 
for all < a < b < oo 

< m((a, b)) < oo. 

For simplicity, it is assumed that m does not have atoms. An impor- 
tant fact is that X has a jointly continuous transition density p(t; x, y) 
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with respect to m, i.e., 

P x (X t eA) = / p(t;x,y)m(dy), 

J A 

where A is a Borel subset of R+. Moreover, p is symmetric in x and 
y, that is, p(t;x,y) = p(t;y,x). The Green or the resolvent kernel of 
X is defined for A > as 

POO 

R\(x,y) = dt e- xt p(t;x,y). 



(ii) Scale function 5 is an increasing and continuous function which can be 

defined via the identity 

where -ff denotes the first hitting time, i.e., 

H y := inf{t : X t = y}, y> 0. 

We normalize by setting 5(0) = 0. Due to the recurrence assumption 
it holds 5(+oo) = +oo. Recall that {S(X t /\H ) '■ t > 0} is a continuous 
local P^-martingale for every x > (see, e.g., Rogers and Williams 
[47| p. 276). It is easily proved that S(X) = {S(X t ) : t > 0} is a 
(recurrent) diffusion taking values in R + . The scale function associated 
with S(X) is the identity mapping x *— > x, x > 0, and we say that S(X) 
is in natural scale. Clearly, also for S(X) the boundary point is 
instantaneously reflecting. Using the Skorokhod reflection equation it 
is seen that S(X) is a P x -submartingale (cf. Meleard |34j Proposition 
1.4 where the semimartingale decomposition is given in case there are 
two reflecting boundaries). 

(iii) The infinitesimal generator of X can be expressed as the generalized 
differential operator 

g = — — 

dm dS 

acting on functions / belonging to the appropriately defined domain 
V{Q) of Q (see ltd and McKean [22], Freedman [12], Borodin and 
Salminen [6]). In particular, since is assumed to be reflecting then 
/ e V(Q) implies that 

x-\o b(x) — 5(0) 
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(iv) The distribution of the first hitting time of a point y > has a P x - 
density: 

P x {H y € dt) = f xy {t)dt. 

This density can be connected with the derivative of a transition den- 
sity of a killed diffusion obtained from X. To explain this, introduce 
the sample paths 

x(y) . = \ x ti l< H v 
1 ' \d, t>H y , 

where d is a point isolated from R + (a "cemetary" point). Then 
{xj; : t > 0} is a diffusion with the same scale and speed as X. Let 
p denote the transition density of X^ with respect to m. Then, e.g., 
for x > y 

fxy{t) = lim —- 2) 

zly S{z) - S(y) 

For a fixed x and y, the mapping t i— > f X y(t) is continuous, as fol- 
lows from the eigen-differential expansions and discussion in Ito and 
McKean p. 153 and 217 (see also Kent [23], [25]). Recall also the 
following formula for the Laplace transform of H y 

E * - tiff < 3 ' 



which leads to 



/•oo 

/ m(dx)E x (e~ aH y) 
Jo 



(v) There exists a jointly continuous family of local times 

{4 y) :t>0,y>0} 
such that X satisfies the occupation time formula 

dsh(X s )= r h(y)L^m(dy), 



where h is a nonnegative measurable function (see, e.g., Rogers and 
Williams [47J 49.1 Theorem p. 289). Consequently, 

1 f l 

L t =1 ] m — 77 T T^vT / l(y-S,y+8)(X s )ds. 

5io m((y - 6,y + 5)) J yy ,a ; 
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For a fixed y introduce the inverse of Z>) via 

:= inf{s : > I}. 

Then = {r { e y) : £ > 0} is an increasing Levy process, in other 
words, a subordinator and its Levy exponent is given by 

B y (exp(-Ar, (?/) )) = exp (-i/R x (y, y)) 

/•oo 

= exp(-W u^(dv)(l-e- Xv )), (4) 
Jo 

where is the Levy measure of . The assumption that the speed 
measure does not have atoms implies that does not have a drift. 
In case y = we write simply L, r and v. 

1.3 Assuming that X is started from we define for t > 

Gt := sup{s < t : A", = 0} and D t := mf{s > t : X s = 0}. (5) 

The last exit decomposition at a fixed time i states that for u < t < v 

P {G t G dn,X 4 G dy, A G du) 

= p(«;0,0) / y0 (i - u) /j/o(« - dudvm(dy). (6) 

In fact, this trivariate distribution is only the skeleton of a more complete 
body of processes: 

{X u :u<G t }, {X Gt+v :v<t-G t }, and {X t+U :v<D t -t} (7) 

the distributions of which we now characterize following Salminen [39J . For 
general approaches; see Getoor and Sharpe [15], [Hj, and Maisonneuve [33] . 

Let x,y G R+ and u > be given. Denote by (X x > u > y , P x ,u,y) the 
diffusion bridge from x to y of length u constructed from X, i.e., the measure 
P x ,u,y governing X x ' u ' v is the conditional measure associated with X started 
from x and conditioned to be at y at time u. The bridge X x ' u,y is a strong 
non-time- homogeneous Markov process defined on the time axis [0, u). For 
the first component in ([7|), we have conditionally on Gt = u 

{X s : < s < G t } = {X° s ' u '° : < s < u} (8) 

For the second component in (|7|) consider the process X^ as introduced 
in (iv) above with y = 0. We write simply X instead of X^°\ For positive 
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x and y let X x ' u ' v denote the bridge from x to y of length u constructed, 
as above, from X. The measure P x , u ,y governing X x ' u ' y can be extended by 
taking (in the weak sense) 

We let X°' u ' y denote the process associated with Po,u,y Then, conditionally 
on Gt = u and X t = y, 

{X Gt +s :0<s<t-G t } = {X°/- U >y :0<s<t-u}. (9) 

For the final part in ([7|), by the Markov property, we have conditionally 
on X t = y 

{X t+S :s<D t -t} = {X s :s> 0}, (10) 

where Xq = y. 



2 Two descriptions of the Ito measure 
2.1 Description due to Ito and McKean 

We discuss the description of the Ito measure n where the excursions are 
studied by conditioning with respect to their lifetimes. Let X be as in section 
1.3 and p(t; x, y) its transition density with respect to the speed measure, in 
other words, 

P x (X t G dy) =P x (X t £dy;t< H ) = p(t; x , y) m{dy) . 

The Levy measure v of r is absolutely continuous with respect to the 
Lebesgue measure, and the density - which we also denote by v - is given by 

v(y) := u(dv)/dv = limlim =: p^v;0,0). (11) 

In Section 1.3 we have introduced the bridge X x,t ' y and the measure 
P x ,t,y associated with it. The family of probability measures {P Xj t, y ■ x > 
0, y > 0} is weakly convergent as y j thus defining P x ,t,o for all x > 0. 
Intuitively, this is the process X conditioned to hit at time t. Moreover, 
letting now x { we obtain a measure which we denote by Po,t,o which 
governs a non-time homogeneous Markov process X '*' starting from 0, 
staying positive on the time interval (0, t) and ending at at time t. 
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Theorem 1. a. The law of the excursion life time £ under the ltd excursion 
measure n is equal to the Levy measure of the subordinator {ti}i>q and is 
given by 

n(( G dv) = v(dv) = p^ (v; 0, 0) dv. (12) 
b. The ltd measure can be represented as the following integral 

POO 

n(de)= / n(Cedt;)Po,„,o(<fe). (13) 
Jo 

Moreover, the finite dimensional distributions of the excursion are charac- 
terized for < t\ < t2 < ■ ■ ■ < t n and x% > 0, i = 1,2, ... ,n by 

n(e tl G dxi,e t2 £ dx 2 , . . . , e tn £ dx n ) 

= m(dxi) f Xl0 (ti) p(t 2 - ti;xi,x 2 )m(dx2) (14) 
x . . . p(t n -t n -x; x n -\, x n ) m{dx n ). 

In particular, the excursion entrance law is given by 

n(e t £ dx) = m(dx) f x o(t), 

and it holds 

n(C>t)= / n(e t €dx)= / m(dx) f x0 (t). (15) 
jo Jo 

Combining the formulas (|12p and (|13p with the last exit decomposition 
leads to a curious relation between the transition densities p and p^ . 

Proposition 2. The functions p(t; 0,0) and p^ (t;0, 0) satisfy the identity 

*t r co 

dup(u;0,0) dv p\v;0,0) = 1. (16) 



'0 Jt-u 

Proof. From (|12p and ()15p we may write 

poo roo 

/ ^p T ( W ;0,0)=n(C>t)= / m(dx)f x0 (t). 
Jt Jo 

Consequently, identity (|16p can be rewritten as 



t poo 

dup(u;0,0) m(dx)f x0 (t- u) = 1, (17) 

but, in view of the last exit decomposition ([6]), identity (|17p states that the 
last exit from when starting from takes place with probability 1 before 
t, in other words, 

Po(G t < t) = I, 

which, of course, is trivially true. □ 



8 



Remark 3. For another approach to Fffl\) notice that it follows from 
and FTl]) 

^ = P„ p r M , 0) (i-e--). 

Hence, from the definition of the Green kernel, 



1= du e- A Xu;0,0) / dv p\v; 0, 0)(1 - e~ Xv ). (18) 
Jo Jo 



OO /"OO /'OO 

du e~ Xu p{u;0,0) / dv e~ Xv ds p ] '(a; 0,0) 

/•OO /'OO 

/ due- A(u+1,) p(u;0,0) / dsp T (s;0,0), 



o 



Consequently, 
1 

A = 



/0 JO Ji) 

/rom which i!6\) is easily deduced. 

2.2 Description due to Williams 

In the approach via the lengths of the excursions the focus is first on the time 
axis. In Williams' description (see Williams [52], and Rogers [45 J, [46]) the 
starting point of the analysis is on the space axis since the basic conditioning 
is with respect to the maximum of an excursion. To formulate the result, 
let for e G E 

M(e) := sup{e t : < t < C(e)}. 

The key element in Williams' description is the diffusion obtained by 
conditioning X not to hit 0. We use the notation and for the measure 
and the expectation associated with X'. To define this process rigorously 
set for a bounded F t G Ct, t > 0, 

El(F t ) := lim B x (F t ;t<H a \H a <H ) 

aT+oo 

lim E x( F f> t < H a<H ) 



aT+oo P x (H a < Ho) 

E x (F t S(X t );t<H a AH ) 

= llm FT~\ ' 

aT+oo b\ x ) 

where the Markov property and formula ([pj for the scale function are ap- 
plied. The monotone convergence theorem yields 
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in other words, the desired conditioning is realized as Doob's /i-transform 
of X by taking h to be the scale function of X. It is easily deduced that the 
transition density and the speed measure associated with X^ are given by 

p\t;x,y) := , m) (dy) := S{yf m{dy). 

We remark that the boundary point is entrance-not-exit for X^ and, there- 
fore, X T can be started from after which it immediately enters (0, oo) and 
never hits 0. 

Theorem 4. a. The law of the excursion maximum M under the ltd ex- 
cursion measure n is given by 

n(M > a) 



S{a\ 

b. The ltd excursion measure n can be represented via 

/•oo 

n(cfe) = / n(M G da) Q^ide), 
Jo 

where Q(*> a ) is the distribution of two independent X^ processes put back to 
back run from until they first hit level a. 

As an illustration, we give the following formula 

n(l-exp(- [ C dsV(e s ))) 
Jo 

= ^°°n(Meda) ^1- (Ej(exp(- j"" duV(ou u ))fj ^ . 
If V > 0, this quantity is the Levy exponent of the subordinator 



{I 



ds V(X S ) : i > } , 

o 



that is, 



E \ exp ( —a J ds V(X S 



= exp < — £ n ( 1 — exp 



-a I dsV{t 
Jo 
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Comparing the descriptions of the ltd excursion measure in Theorem [T] 
(in particular formula (|14p ) and in Theorem H] hints that the processes X 
and have, in addition to conditioning relationship, also a time reversal 
relationship. This is due to Williams [52J, who particularized to the case of 
diffusions the general time reversal result, obtained by Nagasawa [36J. See 
also [33] p. 313, and [6] p. 35. 

Proposition 5. Let for a given x > 

A x := sup{i : u(t) = x} 
denote the last exit time from x. Then 

{X s : < s < H } = {X[ x _ s :0<s< A x }, (19) 
where X$ = x and X = 0. 

3 Stationary excursions; Bismut's description 

Consider the diffusion X with the time parameter t taking values in the 
whole of R. In the case m(R + ) < oo the measure governing X can be 
normalized to be a probability measure. Indeed, in this case the distribution 
of X t is for every t £ R defined to be 

P(Xt E dx) = m(dx)/m(R+) =: fh{dx). 

Recall from ([5]) the definitions of Gt and Dt, and introduce also At := 
D t - G t . 

Theorem 6. Assume thatm(R+) < oo. Then the joint distribution oft — Gt 
and Dt — t is given by 

roo 

P(t - G t G du, D t - t e dv)/dudv = / m(dy) f y o{u) f y o(v) 

Jo 

= v(u + v)/m(R + ). 

Consequently, for At it holds 

P(A t G du)/du = uv{u)/m(R, + ). (20) 

Moreover, the law of the process {Xc t +v ■ v < At} is given by 

C(e)n(de)/m(R + ), (21) 

where n(de) is the ltd measure as introduced in Theorem [I] and [7] and C 
denotes the length of an excursion. 
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Proof. The density of (t — Gt, Dt — t) is derived using the time reversibility 
of the diffusion X, i.e., 



{X t : t G R} = {X_ t : t G R}, 



and the conditional independence given Xt. The fact that the density can 
be expressed via the density of the Levy measure is stated (and proved) in 
Proposition [T2l below, see formulas ([30]) and (j3Tj) . To compute the distribu- 
tion of At is elementary from the joint distribution of t — Gt and Dt — t. For 
these results, we refer also to Kozlova and Salminen [28]. The statement 
concerning the law of {Xc t + V '■ v < At} has been proved in Pitman [37] (see 
Theorem p. 290 point (iii) and the formulation for excursions on p. 293 
and 294) - all that remains for us to do is to find the right normalization 
constant, but this is fairly obvious, e.g., from the density of Aj. □ 

If m(R+) = oo the measure associated with X is still well-defined but 
"only" cr-finite. In this case, the distribution of Xt is plainly taken to be to. 
From (|21h it is seen that we are faced with a representation of the Ito measure 
via stationary excursions valid in both cases m(R+) < oo and m(R+) = oo. 
We focus now on this representation as displayed in (|22p below, and present 
a proof of the representation using the diffusion theory (this provides, of 
course, also a proof of (|2ip ). We remark that in Pitman [37J a more general 
case concerning homogeneous random sets is proved, and, hence, it seems 
worthwhile to give a "direct" proof in the diffusion case. 

Theorem 7. Let F be a measurable non-negative functional defined in the 
excursion space E. Then up to a normalization 



In particular, the process {Xc t + S '■ < s < At} conditionally on At = v 
is identical in law with the excursion bridge X°' v '° as introduced in Section 



Proof. Without loss of generality, we take t = 0. From (|20p we have 



n(F(e)) = E 



( -J- F (X Gt+s : < s < A t ) 



(22) 



2.1. 




Therefore, it is enough (cf. Theorem [T]) to prove that 



n(F(e) | C = u) = E (F {X Go+s : < s < A ) | A = u) . 



(23) 
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P (-G G dv, X G dy | A = u) = fy,o( v ) fy,o(u - v) m(dy)dv. (24) 



Define for < s\ < s 2 < ■ ■ ■ < s n < u 

A 1>n := {X Go+Sl G dyi, . . .,X Go+Sn G dy n }, 

and consider 

poo pu 

E {Ax, n \A Q = u)= / E (Ai,„, -Go G dv, X £dy\A = u) 

J y=0 Jv=0 
poo pu 

= / E (A 1>n I A = u, G = -w, X = y) 

Jy=0 Jv=0 

xP (-G G dv, X G dy \ A = n) . 

From the description of the process X, the conditional independence and 
the equality of the laws of the past and future given Xq, and using formula 
ZU1) we obtain 

1 

uv(u) 

Letting k be such that — v + < t < —v + Sk+i, if any, we write applying 
again the conditional independence 

E (Ai }n I A = u, G = -w, Xq = y) 

= E (Ai ;k A k+1>n | A = u, G = -u, X = y) 

= E (Ai )fc | G = -v, X = y) E (A fc+ i !n | D = u - v, X = y) . 

Recall from Introduction section 1.2 (iv) the notation X for the diffusion X 
killed when it hits 0. As in section 1.3 of Introduction we may construct the 
bridge X y,v, ° starting from y having the length v and ending at 0. We let 
Py,i),o denote the measure associated with X y ' v '°. With these new notations, 

E (A hk | G = -v, X = y) 

= Py.u.o {u v -s k G dy k , ... , io v - Sl G dyi) 

p(v - s k ;y,yk)m{dy k )p{s k - s k -x;y k ,y k -i)m(dy k -x) 



fyo{v) 

x ... p(s 2 - s 1 ;y 2 ,yi)m(dyi) f yi o( s i) 

and 

E (A k+1)n | D = u - v, X = y) 

= Py,v,o {u Sh+1 - v G dy k+ i, . . . ,uj Su - v G dy n ) 
■p( s k+i-v;y,y k+ i)m(dy k+1 ) 



fyo(u - v) 

x p(s k+ 2 ~ Sfc+i; Vk+i, Vk+2) m(dy k+2 ) 

x ... p(s n - s n _i; y n -i,y n ) m(dy n ) f y n o(u s r 
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Using now (|24p and formulas above we have after some rearranging and 
applying the symmetry of the transition density p 

E (A ltn | A = u) 
1 

= TT m ( d y±) fyio(si)p{s 2 - si; yi,y 2 ) m(dy 2 ) 

uv[u) 

pu poo 

x... dv m(dy)p(v - s k ;y k ,y)p(s k+1 -v;y,y k+1 ) 
Jo Jo 

x . . . p(s n - s n -i;y n -i,y n ) m(dy n ) f Vn 
Performimg the integration yields 

E (A 1>n | A = u) 

= m(dyi) f yi o{si)p{s 2 - sf, y t , y 2 ) m{dy 2 ) 

and this means that (|23|) holds completing the proof. □ 

Remark 8. The formula A22\) was derived for Brownian motion by Bis- 
mut J3j/. The connection with the Palm measure and stationary processes is 
discussed in Pitman [31]. In fact, Bismut describes in the Brownian case 
the law of the process {Xc t+S : < s < At} in terms of two independent 
3-dimensional Bessel processes started from and killed at the last exit time 
from an independent level distributed according to the Lebesgue measure (see 
JBS and for details). 

4 On the excursion straddling an independent ex- 
ponential time 

In the literature one can find several papers devoted to the properties of 
excursions straddling a fixed time t; first of all, Levy's fundamental paper 
|31j . which contains a lot about the zero set of Brownian motion, its (inverse) 
local time, excursions, and so on. See also Chung [7] starting from Levy's 
paper [31], Durret and Iglehart [9], and Getoor and Sharpe [16], [TTJ. In fact, 
the last exit decomposition ([6|) lies in the heart of these studies (see Getoor 
and Sharpe [E], [Il|)- However, it seems to us that excursions straddling 
an exponential time are not so much analyzed. Here we make some remarks 
on this subject. 
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Let T be an exponentially distributed random variable with parameter 
a > 0, independent of X, and define 

G T := sup{s <T :X S = 0}, D T : = inf{s > T : X s = 0}, 

and 

A7 1 := — Gt ■ 

The Levy exponent of the inverse local time at is denoted by <£, in other 
words, 

E (exp(-Ar £ )) = exp {-I $(A)) 
Recall the relation (cf. Q with y = 0) 

$(A) J R A (0,0) = 1. (25) 
4.1 Last exit decomposition at T 

We begin by discussing the last exit decomposition at the exponential time 
T. 

Theorem 9. (i) The processes 

{X u :u< G T } and {X Gt+v : v < A T } 

are independent. 

(ii) The law of {X u : u < Gt} may be described as follows: 

(a) Lt ■= Tq t is exponentially distributed with mean l/$(a), 

(b) The process {X u : u < Gt} conditionally on Lt = £ is distributed as 
{X u : u < ti} under the probability 

exp (-are + £<&(a)) Pq- 



(iii) The law of the process {Xg t+v : v < A^} is given by 

^(l-e-^))n(d £ ). (26) 

where n(cfe) is the ltd measure associated with the excursions away from 
for X and £ denotes the length of an excursion. 
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Proof. Let F\ and F2 be two nonnegative functionals of continuous functions 
and consider 

E (F 1 (X U :u<G T ) F 2 (X Gt+v : v < A T )) 

POO 

= a dt e" at E (F^Xv : u < G t ) F 2 (X Gt+v : v < A t )) 
Jo 

= qE ( f'dt z~ at Fi{X u : u < n-) F 2 (X Te _ +v : v < r t - r<_) 

= E (Y die-^F^Xu :u<n)\ 

x /" n(ck) (l - e- a « e )) F 2 (e s : s < C(e)) 

= *(a)E Qf dl<r aT *F x {X u :u<n)\ 

r /i_ e -«C(e)\ 
x y n(cte) I J F 2 (e s : a < C(e)). 

where the third equality is based on the properties of the Poisson random 
measure associated with the excursions (see Revuz and Yor |43| Master 
Formula p. 475). □ 



Remark 10. Notice that letting a — > in \26\) and using 

lim = lim a R a (0, 0) = 1/WR+) (27) 

a-^0 <P(a) a^O 

yield the probability law of the excursion straddling a fixed time in the sta- 
tionary setting, cf. h21\) in Theorem® 

As a corollary of Theorem[9]we have the following results which show that 
after conditioning the quantities do not depend on a, and in this context a is 
entirely "contained" in Gt and At- The formulas should be compared with 
©) ©, and (fTUj) . The distributions of Gt and At are given, respectively, 
in (HOT) and (EH) below. 



Corollary 11. For any nonnegative functionals F\ and F 2 of continuous 
functions it holds 

E {F X {X U :u<G T )\G T =g) = ^o {Fi(X u :u<g)\X g = 0) 

= E ,g,o{F 1 (u Ju :u<g)) (28) 
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and 

E (F 2 (X Gt+v : v < A T ) \ A T = h) = E (F 2 (X v : v < h) \ X h = 0) . 

= E 0A0 (F 2 {uj s :s<h)) (29) 

Proof. The statement (|28p can be obtained from the corresponding result for 
fixed time as presented in (JS]). Also (|29|) can be derived from the fixed time 
result but we prefer to present here a proof based on the Master Formula. 
For this consider for 5 > 

E (F 2 (X Gt+v : v < A T ) 1{/i<a t </i+5}) 

r /i_ e -«C(s)\ 
= y n (&) I J *2(e« : a < COO) 1 {h<((e)<h+8}- 

Using the description (|13p of the ltd excursion law we obtain 

E (F 2 (X Gt+v : v < A T ) l{h<A T <h+5}) 
r-h+S (y_ e -au\ _ 

duv(u) Eo,«,o(^(w s : s < u)). 

h V / 

Applying the explicit form of the distribution of At as given in (|35p and 
letting 5 j leads to {22} • □ 

4.2 On the distribution of (G T , D T ) 

In this section the distributions of T — Gt, Dt — T and At := Dt — Gt are 
studied in detail. 

Proposition 12. The joint distribution ofT — Gt and Dt — T is given by 
P (T -G T edu,D T -T e du) 

/•oo 

= dudvaR a (0,0)e- au / m(dy) f y0 (u) f y0 (v). (30) 

JO 

= ae-»M, + ,) ^ 
$(a) 

In particular, 

poo 

v(u + v)= m(dy) f y0 (u) f y0 (v). (32) 

JO 
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Proof. From ([6]), 

P (* - G t G du, D t - t £ dv) 

= dudvp(t - u;0,0) l{ u <t,v>0} / m(dy) f y0 (u) f y0 (v), 

and, hence, 

P (T — Gt € du, Dt — T £ dv) 

= dudva / dte~ at p(t - it; 0,0) / m(dy) f y0 (u) f y o(y) , 



o 





from which (|3Q|) follows. To derive (|3ip. we apply again the Master Formula 
(see Revuz and Yor [33] p. 475). For this, let (u,v) t— > (p(u,v) be a non- 
negative and Borel measurable function and define 

Q(y) :=V^(T -G T ,D T -T)). 

Letting ti denote the inverse of the local time I at we have 

Q(tp) =E J] ip(T - Ti^,Ti - T) l{ n _ <T <r e } 

= Eo (j^ <p(T -Tt,z + rt- T) l{ n<T<n +z} v ( z ) dzd ^j , 

since {((,Tg) : I > 0} is a Poisson point process with Levy measure didu, 
and T is independent of \ji : £ > 0}. Apply next that T is exponentially 
distributed to obtain 



Q(p) = E UT_ 



r-Tf + z 

v{z)dzdl I dtae~ at f(t — re,z + re — t) 



} e -a(x+r e ) ^ z _^ ± dxdzM ] ^ 



'K. 2 -/r, 

= qE I / i/f«) 

where we have substituted x = t — T£. Furthermore, setting y = z — x yields 
Q(<p) = aE (^J s + x) e - Q(a;+T * ) ^(x, y) cfedyd^ 

= aE e " 073! ^ / 2 ^(x,y)e- ax u(y + x)dxdy, 

and (|3ip follows now easily from (UJ. The equality (|32p is an immediate 
consequence of (f5U|) and (|3T|) . □ 
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Corollary 13. 1. The densities for T — Gt, Dt — T, and At are given, 
respectively, by 

P (T — Gt G du)/du = e~ au [ u(z)dz, (33) 

Ju 

P (D T -Te = -^-e av r e- az u(z)dz, (34) 



and 



P (A T G da)/da = (1 l™^ ' ^ 

9(a) 



2. The joint density ofT — Gt and At is 

01 

Po(T-GTedu,A T eda)/duda = —-e- au u(a), u < a. (36) 

<P(a) 

3. The density ofT — Gt conditionally on At = a is 

OL 

P (T-G T £du\AT = a)/du = e~ au , u < a. (37) 

1 — e~ aa 

Proposition 14. The joint Laplace transform of Gt and Dt is given by 

t? ( r n $(72 + «) ~ $(72) f ^ 

E (exp (-71 G T - 72-Dr)) = — : ; — v— • (38) 

<P(7l + 72 + a) 



In particular, 



Eo = »(^«)-«M (39) 

* a 



and t/ie random variables Gt and At are independent. The density of Gt 
is given by 

P (G T G du)/du = $(a)e- Q >(u;0,0). (40) 

Proof. The formula (|40p for the density of Gt is obtained from (|6|) by in- 
tegrating. The independence of Gt and At follows immediately from 
To derive the joint Laplace transform of Gt and Dt, consider 

E (exp (-71 G T - 72 -Dt)) 

= y dtae~ at J e' JlU ' j2V P (Gt G du, A G dv). 
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Applying the last exit decomposition formula © yields 
E (exp (-71 G T - 72 D T )) 

f'QO f-V 

due-^ u p{u; 0,0) / dve~^ v / dtae~ at 



x / m(dy) fyo(t - u) fyo(v - t) 
Jo 

poo poo pa+u 

/ due-^ u p(u; 0,0) / da e -^ {a+u) / dt a e~ at p{u- 0, 0) 

/>oo 

x / m(dy) f y0 (t - u) f y0 (a + u - t) 
Jo 

due-^ 1+ ^ )u p(u; 0,0) / dae^ / dba e - a(fe+M) p(n; 0, 0) 
Jo Jo 

poo 

x / m(dy) f y0 (b) f y0 (a - b) 
Jo 

X) /*oo /*oo 

a / rf ne - (71+72+a) >(u;0,0) / m(dy) / dae" 72a 



x / dbe- ab f y0 (b)f y0 (a-b) 
Jo 

= a« 7l+72+a (0, 0) m(dy) E y (e^ )* ) E y (e" 72 " ) . 
To proceed, we have 

J™m{dy)E y (V (72+a) ^°) By (e~^ H °) 

m(dy)Ry 2+a (y,0)Ry 2 (y,0). 



i? 72+Q (0, 0)i? 72 (0, 0) J 
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The integral term in this expression can be evaluated: 
m(dy) R a+12 (y, 0)i? 72 (y, 0) 

) f'OO P oo 

m{dy) I dte~ {a+ ^ 2)t p(t;y,0) / ds e" 72S p(s; y, 0) 
Jo Jo 

f'OO poo 

/ dte~ {a+r2)t dse-^ s p{t + s; 0,0) 
Jo io 

/'OO 

dte -(«+7 2 )t / due-^-Qpiu; 0,0) 



oo i — an 

dne" 72n p(«;0,0), 

o « 

= i(i? 72 (0,0) -i? 72+Q (0,0)). 

where the Chapman-Kolmogorov equation and the symmetry of the transi- 
tion density p is applied, and by ([25]) this completes the proof. □ 

Remark 15. 1. From Proposition it is seen that the density of At can 
also be written in the form 

a i , , , ; „ _„,, 



P (A T G da)/da = — / m{dy) / dbe~ ab f y0 (b) f y0 (a - b), 



which taking into account 135\) leads to the identity 
(1-e- 



v(a)= / m(dy) / dbe- Qb f y0 (b) f y0 (a - b). 
a Jo Jo 

Let here a —* to obtain 

f°° f a db 

u(a)= m{dy) - f y0 (b) f y0 (a - b) . (41) 
jo jo a 

It is interesting to compare this expression with the following one obtained 
from (SB) 



u{a) = / m(dy)f yQ (b)f y0 (a-b). (42) 
J o 

The fact that the right hand sides of |^1| ) and fl^i?| ) do no£ depend on b can 
also be explained via the Chapman-Kolmogorov equation. 
2. We may study distributions associated with Gt, Dt and At in the station- 
ary case, i.e., if m(R+) < oo, by letting a — * 0, as observed in Remark 1 lffl 
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From Proposition UM o,nd Corollary [75] we deduce the following results take 
t = 0: 

P(— Go G du, Dq G dv)/dudv = — — - — - v(u + v). 

m(R+) 

i r 00 

P(-G G du)/du = P(A) G du)/du = — / dv, 

m(bl + ) J u 

P(Ao G da) Ida = — — - — - a via). 

m(R+) 

Moreover, letting Zt '■= (T — Gt)/At then (Zy, Ay) converges in distribu- 
tion as a — > to ([/, A), where U and A are independent with U uniformly 
distributed on (0,1) and A is distributed as Aq fc/. T/ieorem[^j. 



4.3 Infinite divisibility 

In the paper by Bertoin et al. [2] it is proved that the distribution of Ay 
for a Bessel process with dimension d = 2(1 — a),0 < a < 1, is infinitely 
divisible (in fact, self-decomposable) and the Levy measure associated with 
this distribution is computed. In this section we show that the distribution 
of At is infinitely divisible in general, i.e., for all regular and recurrent 
diffusions. Moreover, we also prove that the distributions oIT—Gt and D?— 
T have this property. The key to these results is the Krein representation 
of the density of the Levy measure v (see Knight [26], Kent [25], Kiichler 
and Salminen [30], and, in general on Krein's theory of strings, Kotani and 
Watanabe [27], Dym and McKean [10] ) according to which 

/■oo 

v{a) = / e~ az M(dz), (43) 
J o 

where the measure M has the properties 

f°° Midz) , f°° M(dz) 

/ — - — '- < oo and / — - — - = oo. 

Jo z ( z + !) Jo z 

Theorem 16. The distributions of T — Gt, Dt — T and At are infinitely 
divisible. 

Proof. As seen from (j33j) . ()34|) . and (|35p . the intrinsic term in the densities 
of T — Gt, Dt — T and A^ is the density v(a) of the Levy measure of 
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the inverse local time at 0. We consider first the distribution of T — Gt- 
Applying the Krein representation (|43p in (|33|) yields 

P (T — Gt G du)/du = a R a (0, 0) e~ QU da M(dz) e~ az 

Ju JO 



« , e -au r M(dz) ^ uz 



a r M(dz) c _ (a+z)u 



$(a) 

(a + z)e- (Q+2) "M a (dz), 



with 

/ , \ a M(dz) , 1A . 
M a {dz) = — \> 44 
£>(a) z(a + 

The claim of the theorem follows now from the fact that mixtures of expo- 
nential distributions are infinitely divisible, (see Bondesson [5]). For Dt — T 
we compute similarly from (|34p via the Krein representation 

P (D T -Te dv)/dv = e av f°° Q- aa v{a)da. 

^\OL) j v 

ze- zv M a (dz). 



To analyze the distribution of At we use the Krein representation in fl35f) 
to obtain 

1 f°° / \ 
P (A T e da) Ida = — — / ( e~ za - g — ( Q + Z ) a j M(dz). (45) 

n«) Jo v y 

Notice that for a > 

/( a; z a ) = Z ( a + Z ) ( e -za _ e -(a+z)a 

a \ 

is a probability density as a function of a. In fact, letting T± and T2 be 
two independent exponentially distributed random variables, with respective 
parameters z and a + z, then the sum T\ + T2 has the density f(a; z, a). In 
particular, the distribution of T% + T2 is a gamma convolution (which, by 
definition, is the law of finite sum of independent gamma variables). Next 
we notice that letting 

z(a + z) _ 2 

Tl z a as := x dx, z < x < a + z 

a 
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we may represent the distribution of T\ + T2 as a mixture of Gammas- 
distributions as follows 

POD 

f(a;z,a)= x 2 ae~ xa IL z>a (dx). (46) 
Jo 

Combining the representation (|46p with (|45p yields 

P (A T e da)/cfo = ^-J'°° ^4 
$(a) Jo z(a + z) 

/>oo 

= / x 2 ae-"n Q (ds), (47) 
Jo 

where II Q is a probability measure on R + given for any Borel set A in R + 
by 

poo 

R a (A) = / M a (dz)U z>a (A). (48) 
Jo 

The claim that the distribution of A-p is infinitely divisible follows now 
from (|47l) by evoking the result that mixtures of Gamma(2)-distributions 
are infinitely divisible (see Kristiansen [29] ). □ 

Remark 17. 1. Recall from Bondesson that a probability distribution 
F on R+ is called a generalized gamma convolution (GGC) if its Laplace 
transform can be written as 

J™ e- sa F(da) = exp (-ps + log (j^A U(dt)\ , (49) 

where jjL > and U is a measure on (0, 00) satisfying 

I \\ogt\U(dt) < 00 and f U (dt) ^ ^ 
7(0,1] J(i, 00) * 

is known see f^j Theorem 4-1-1 p. 49 that if (3 is the total mass ofU then 
the distribution F in fi49\ ) is a mixture of Gamma ((3) -distributions 
2. The distribution of the length At of an excursion straddling a fixed time 
t for a stationary diffusion (with stationary probability distribution) is given 
in Theorem^ \20\) as 
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Also in this case the distribution of At is a mixture of Gamma(2)-distributions 
and, hence, it is infinitely divisible. In fact, 

roo 

P(A t Gda)/da= z 2 ae~ za M(dz). 
Jo 

where the probability measure M is given in terms of the Krein measure M 
via 

M(dz) = M(dz)/(m(H + )z 2 ). 

5 Case study: Ornstein-Uhlenbeck processes 

In this section we give some explicit formulas for excursions from to 
associated with Ornstein-Uhlenbeck processes. It is possible to obtain such 
formulas due to the symmetry of the Ornstein-Uhlenbeck process around 
0. Analogous results for excursions from an arbitrary point x to x are less 
tractable. 

5.1 Basics 

Let U denote the Ornstein-Uhlenbeck diffusion with parameter 7 > 0, i.e., 
U is the solution of the SDE 

dllt = dB t — jUtdt with Uq = u, 

and most of the time, but not always, we take u = 0. Recall that the speed 
measure and the scale function of U can be taken to be 

2 f x 2 

m(dx) :=2e~" ,x dx and S(x) := / e lv dy, 

Jo 

respectively. Moreover, see [6] p. 137, the Green kernel of Ornstein-Uhlenbeck 
process with respect to the speed measure is given for x > y by 

x exp f^j D_ xh (-3/^27) . 
where D denotes the parabolic cylinder function. In particular, since 
£U/ 7 (0) = (2 A/(27) r ((A + 7 )/(2 7 ))) 1 , 
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we have, after some manipulations, 



i? A (0,0) = ^ {X J 1] (2^7) r((A + 7)/(27)) 



2^7 

Consequently, using the formula 



r(x) = ^r((x + i)/2)r(x/2) 



we obtain 

r> (r\ r\\ 1 r(A/(2 7 )) 

^ (0 ' 0) = ¥(A) = 4r((A + 7 )/(2 7 ))- (50) 

We remind also that U can be represented as the deterministic time 
change (Doob's transformation) of Brownian motion via 

U t = e" 7i (u + (3 at ), 

where j3 is a standard Brownian motion and at := (e 27 * — l)/2 7 (see Doob 
0). 



5.2 Killed Ornstein-Uhlenbeck processes 

We consider now the Ornstein-Uhlenbeck process killed at the first hitting 
time of 0, and denote this process by U. Let Y be the diffusion on R + 
satisfying the SDE 

dY t = dB t + (y f - 7*t) dt, Y = y>0. 

Recall that Y may be described as the radial part of the three-dimensional 
Ornstein-Uhlenbeck process. In [6] p. 138 the basic properties of such 
processes are presented. In particular, we record that is an entrance-not- 
exit boundary point and the process is positively recurrent its stationary 
distribution being the Maxwell distribution, i.e., the distribution with the 
density proportional to the speed measure of Y, that is, 

m Y (dx) := 2x 2 e" 7 ^ 2 dx, x > 0. 

We remark that there is a misprint in [6] p. 139; the stationary distribu- 
tion in the general case is not a x 2 -distribution but a generalization of the 
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Maxwell distribution. The transition density of Y with respect to its speed 
measure m Y is 

Y( t . _ y^e 3 ^ 2 / 7e -7*(^ + y2) x / lxy \ 

P V27rsinh( 7 t)xy ^ V 2sinh( 7 t) y Sm ^sinh( 7 i)y' 

and can be computed from the transition density of a Bessel process using 
Doob's transform (for an approach via inverting the Laplace transform see 
Giorno et al. [IS]). In Salminen [38] it is proved that 

Px(Ut e dy) = P x (U t edy,t< H ) 

= e-^p Y (t-x,y)^\m Y (dy), (51) 

where ip~/(x) = 1/x is the unique (up to multiplicative constants) decreasing 
positive solution of the ODE associated with Y killed at rate 7: 



2 

From (1511) we obtain 



-u"(x) + ( 7 x^ u{x) = ju(x). 



Proposition 18. The transition density (with respect to its speed measure 
m) of the Ornstein- Uhlenbeck killed at the first hitting time of is given by 

^ e7,/2 exD (^*-+y 2 y\ smh ( 



m x ' y) = jssm exp r J smh v^hw J - <52) 



Combining the expression of the transition density in (|52p with formula 
([2]) yields the distribution of i^o (see also Sato [50] and Going-Jaeschke and 
Yor EE]). 



Proposition 19. The density of the first hitting time of for the Ornstein- 
Uhlenbeck process {Ut} is given by 

fx0it) = v^(sinh( 7 t))3/ 2 6XP r^nhRj • (53) 

5.3 Levy measure of inverse local time and densities of At, 

T — Gt and D T -T 

The density of the Levy measure of the inverse local time at is obtained by 
applying formula (jlip (see also Hawkes and Truman [20]). Moreover, using 
(|50p in formula ^ leads to an explicit expression for the Bernstein function 
associated with the inverse local time at 0. 
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Proposition 20. The density of the Levy measure of the inverse local time 
at is 



7 3 / 2 e^/ 2 (2 7 ) 3 / 2 



.27* 



f (*) = -f^ ~ r = r— ~ r- (54) 

W ^(sinh( 7 t)) 3 /2 v^( e 27t- 1)3/2 ^ ' 



Lei {r£ : £ > 0} he the inverse local time at 0. Then 

E (exp(-A T ,))^,p(-, 4r ^y^» ). 

Next we display the distributions of Ay, T — Gt, and Dt — T. Recall that 
these distributions are infinitely divisible and the densities are expressable 
via the density of the Levy measure, as stated in Corollary [13] formulae (|33p 
and (|34p . and in Theorem 1161 To simply the notation, we take 7 = 1. 

Proposition 21. Wist/i <I>(a) as in &5(J\) , the distributions of At, T — Gt 
and Dt — T are given, respectively, by 

P (A T G da)/da = - . -= e 2a (e 2a - l) 3/2 , (55) 
<P(a) V 71 " 

Po(r _ Gt e da)/da = ( e 2« _ !) V2 (56) 

and 

P (D T -Te da) /da = ^f- / due - au —e 2u (e 2u - l) 61 . (57) 

*(a) Ja 

5.4 The Krein measure 

As seen in Section 4.3, the Krein representation plays a central role in the 
proof of infinite divisibility of the distributions of T — Gt-, Dt — T, and At- 
Therefore, it seems motivated to compute the measure M (cf. ()43[) ) in this 
representation for Ornstein-Uhlenbeck processes. 

To start with, we give the spectral representation of the transition density 
of p of the Ornstein-Uhlenbeck process killed at the first hitting time of 0. 
Instead of computing from scratch, we exploit the spectral representation 
for p Y (with 7 = 1) as presented in Karlin and Taylor [23] p. 333: 

00 

p Y (t;x,y) =5X1/* e-**LW*)tf)LMtf), (58) 

n=0 
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where {Ln : n = 0, 1, 2, . . . } is the famly of Laguerre polynomials with 
parameter 1/2 normalized via 

[ (4 1/2 V)) 2 m Y (dx) = ^( U+ J)=-- -n,i/ 2 - (59) 

Notice that we consider the symmetric density with respect to the speed 
measure m Y . From (I5ip and ([58]) the spectral representation of p is now 
obtained immediately and is given by 

oo 

p(t;x,y) = E<i/2 e-( 2 " +1 )*x4 1/2) (^)y4 1/2) (y 2 )- (60) 



n=0 



The normalization (|59p coincides with the normalization in Erdelyi et al. [TT] 
(see formula (2) p. 188 where the notation for the norm is h n ). Therefore, 
from [TT] formula (13) p. 189 we have 



n + 1 



n 



4 1/2) (o)= T 5 (6i) 



and, consequently (cf. (|53p ). we obtain the spectral representation for the 
density of the first hitting time of 

oo 

/*>(*) = 5X1/2 e-( 2 " +1 ^4 1/2) (x 2 )4 1/2) (o). 

n=0 

„ oo 

= J-Y e-V^xL^ix 2 ). (62) 

To find the spectral representation for the density of the Levy measure we 
apply formula (|54|) which yields 

"M=4£rt 5 V (2 " +i ". ess) 



7r " — ' \ n 

n=0 



In view of (|43p . we have 



Proposition 22. TTie measure M in the Krein representation of v for the 
Ornstein- Uhlenbeck process is given by 



n=0 

where 5 is the Dirac measure. 



9 00 / I 1\ 
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Notice that 

„( t ) = _? =e -t (!_ e -2t)-3/2 



and, hence, (|63p is obtained also from the MacLaurin expansion of x 
(1 — x) -3 / 2 evaluated at x = e~ 2t . 



Acknowledgement. We thank Lennart Bondesson for co-operation con- 
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